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Abstract 

The process e + e _ —► KiKg has been studied with the CMD-2 detector using about 
950 events detected in the center-of-mass energy range from 1.05 to 1.38 GeV. 
The cross section exceeds the expectation based on the contributions of the p( 770), 
a;(782) and 0(1020) mesons only. 
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1 Introduction 


The investigation of the process e + e” —»• KgKg is important for a number of 
physics problems. Since both 1=0 and 1=1 vector mesons can decay into a kaon 
pair, one can search for excitations of the p(770), cu(782) and 0(1020) by mea¬ 
suring the cross section of the process in the energy range above the 0(1020) 
meson [1], The isovector part of the cross section of the process e + e“ —»• KK 
(both K + K~ and Kg Kg final states should be taken) can be related to the 
t~ —> K~K°u t decay by using conservation of vector current [2], Assuming 
the hypothesis of factorization it can be also used to account for the produc¬ 
tion of kaon pairs in B~ —> D°K~K° decays [3]. Finally, the process under 
study contributes to the total hadronic cross section of e + e _ annihilation, so 
that the values of its cross section are used in the calculation of the hadronic 
contribution to the muon anomalous magnetic moment [4]. In view of the in¬ 
creasing experimental accuracy in the measurement of this quantity [5], any 
significant contribution like that from the process e + e~ —> KgKg should be 
measured with adequate precision. 

Earlier measurements of the cross section performed by the DM1 collaboration 
in Orsay [1] and at OLYA and CMD detectors in Novosibirsk [6,7] were based 
on small data samples and had a systematic accuracy of about 20% or worse. 
Significant progress in the study of the process e + e~ —> KgKg was achieved 
by the SND collaboration [8] at the VEPP-2M collider [9]. The experiment 
was based on integrated luminosity of 6.3 pb _1 . The systematic error of the 
cross section was estimated to be 10% around 1.04 GeV increasing to about 
15% at 1.38 GeV. 

In this work we report on the measurement of the e + e“ — > KgKg cross section 
based on 5.7 pb -1 of data collected with the CMD-2 Detector [10] at the 
VEPP-2M collider from 1.05 to 1.38 GeV. The systematic uncertainty on the 
cross section was about 5% below 1.09 GeV and increased to 10% at 1.38 GeV. 


2 Detector and experiment 


The CMD-2 detector has been described in detail elsewhere [10]. Its track¬ 
ing system consists of the cylindrical drift chamber (DC) [11] surrounding the 
interaction point, and proportional Z-chamber (ZC) [12] for precise polar an¬ 
gle measurement, both also used for the trigger. Both chambers are inside a 
thin (0.38 A 0 ) superconducting solenoid [13] with a field of 1 T. The barrel 
electromagnetic calorimeter [14] is placed outside the solenoid and consists 
of 892 Csl crystals. The muon-range system [15] of the detector, also located 
outside the solenoid, is based on streamer tubes. The endcap electromagnetic 
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calorimeter [16] based on BGO crystals makes the detector almost hermetic 
for photons. 

The data sample used in the analysis was collected in two scans of the center- 
of-mass energy range 1.05 - 1.38 GeV. In the first scan the beam energy was 
increased from 530 MeV to 690 MeV with a 10 MeV step, while in the second 
one it was decreased from 685 MeV to 525 MeV with the same energy step. 


3 Data analysis 


Events of the process e + e —■> K®Kg were detected by using a subsequent 
decay of the Kg meson into a pair of charged pions. 

The following selection criteria were used: 

• There are two oppositely charged tracks from the vertex closest to the beam. 
Track momenta, assuming that tracks are pions, satisfy the conditions: 

pmin _ 20 < P 12 < P" MX + 20, 

where P™ m , P™ ax are minimum and maximum kinematically possible mo¬ 
menta of a pion from the Kg —> decay in MeV/c. 

• The track polar angles are: 


0.95 < 0 1)2 < vr-0.95. 


• The maximum of the track ionization losses is 

,,dE. ,dE ,dE 

ro< “ (( * )l ' ( * )2)<2 ' 2( * w ’ 

where is the ionization loss of a minimum ionizing particle (see 

Fig. la for events from the energy range 1.05 < •>/s < 1.09 GeV). The cut is 
used to remove K + K~ events as well as events of “beam-wall” interactions. 

• The acollinearity angle for pions in the R -93 plane is: 


0.2 < |7T - \ipi — v? 2 11 < 3.0, 


and the space angle between pion tracks is 

^ > 0.5. 

By these criteria events with collinear particles and particles going close to 
each other were rejected. 
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Fig. 1. Distributions used for the selection of e + e~ —> events: a. Particle 

ionization losses versus the invariant mass of two tracks; b. Z-coordinate of the 
vertex versus the distance from the beam point to the vertex in the R — p plane; 
c. Missing momentum (the histogram is simulation, the points are experimental 
events); d. Invariant mass of two tracks. 

• The Z-coordinate of the vertex is less than 7.0 cm and the radius of the 
vertex in the R -92 plane ( Rvrtx ) satisfies the condition: 

0.07 < R V rtx < 1-3 cm, 

as shown in Fig. lb for the energy range 1.05 < a/s < 1.09 GeV. 

In the energy region above the 0(1020), production of a neutral kaon pair is 
often accompanied by emission of a hard photon by initial electrons (“return 
to resonance” effect). The distribution of the missing momentum of two tracks 
defined as P m ,i S = I-P 1 +-P 2 I is shown in Fig. lc for ^fs = 1.07 GeV. In agreement 
with the Monte Carlo simulation, the left peak in Fig. lc corresponds to 
“return to resonance” events, while the right one describes events without 
hard photon emission. In the present work events with “return to resonance” 
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were excluded from analysis by the requirement: 


Prr 


~ '/ E beam ~ m K° 


< 40 MeV/c, 


where 40 MeV/c corresponds to five standard deviations of the experimental 
resolution for the P mis value. The selection criteria are shown in Fig. lc by 
the vertical lines. 


The number of events was determined from a fit of the invariant mass distribu¬ 
tion with a sum of two Gaussians describing the signal and a smooth function 
describing background. For the fitting procedure, the data sample was sub¬ 
divided into three energy bins: 1.05-1.19, 1.20-1.29 and 1.30-1.38 GeV. The 
approximation above was performed in each of the three bins and the param¬ 
eters obtained were later fixed during the approximation at each energy point 
within the corresponding bin. An example of such an approximation in the 
energy range 1.05-1.19 GeV is shown in Fig. Id. 

After background subtraction and application of the cuts described above, 
948±33 KiKg events were selected. 

At each energy point the cross section is determined from the following for¬ 
mula: 

N 

Le{ 1 + S ra d) 

where N is the number of selected events, £ is the detection efficiency, L is 
the integrated luminosity determined from events of Bhabha scattering at 
large angles [17], and (1 + S ra d ) is a radiative correction due to initial state 
radiation [18]. The detection efficiency £ = £ r ec£tri g £geom, where £ rec is the 
reconstruction efficiency, £ trig is the trigger efficiency, and £ ge om is the accep¬ 
tance. 

The reconstruction efficiency (£ rec ) was determined from the experimental 
data [19]. To this end “test” events, in which a K® produced a cluster in the 
Csl calorimeter, were selected. Using the angles of this cluster as well as the 
angles of the clusters produced by pions from the K° s —> 7r + 7r - decay and 
requiring one track in DC, one obtains a clean sample of events for efficiency 
determination. Analysis shows that the reconstruction efficiency is energy in¬ 
dependent. A similar procedure was used to determine the trigger efficiency 
(Et.rig) [19]. The acceptance ( £ ge0 m ) was determined from Monte Carlo simula¬ 
tion. 


The beam energy at each point was evaluated from the value of the magnetic 
held in the dipole magnets [20]. The systematic uncertainty was estimated to 
be ^ = 4 • 10 -4 from the analysis of the long-term stability of energy [17]. 
The number of events, integrated luminosity, detection efficiency, radiative 
correction and cross section at each energy point are listed in Table 1. It is 
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this cross section (the “dressed” one) that should be used in the approxima¬ 
tion of the energy dependence with resonances. For applications to various 
dispersion integrals such as the leading order hadronic contribution to the 
muon anomalous magnetic moment, one should use the “bare” cross section. 
Following the procedure in Ref. [21], the latter is obtained from the “dressed” 
one by multiplying it by the vacuum polarization factor |1 — fl(s)| 2 , where 
fl(s) is the photon polarization operator calculated taking into account the 
effects of both leptonic and hadronic vacuum polarization. 

Figure 2 shows the energy dependence of the cross section obtained in this 
work together with the results of our study of the 0 —► K ( [Kg decay from [22] 
(the data in the 1.00 - 1.04 GeV energy range). Also shown are the results of 
the previous experiments that studied the process e + e^ —> K^Kg above the 
0 meson region. Good agreement between the results of all measurements is 
observed. 

The energy dependence of the cross section was approximated in the frame 
of Vector Dominance Model [23] with the contributions of the p( 770), a;(782) 
and 0(1020) mesons according to the following SU(3) based formula: 


u(s) 


1 q 3 (s) r ^\/u(0 -> K 0 L K° s )m^ 

s 5/2 q 3 (ml) D^s) 



-> e + e~)B((j) - 

- KIK%) 

D w (s) 


y/Y<t,Y p mlm 3 p 6iTB(p - 

+ e + e~)B((p ~ 

* KK° S ) 

D„M 



where cr(0 —> K°Kg) = 127rR(0 —■> e + e )R(0 —> KgKg)/m^ is the only free 
parameter, determined from the fit, q(s) = \Js /4 — m 2 K o is the neutral kaon 
momentum, and D v (s ) = — s + 7 v /iFy(s). The energy dependence of the 

total width for a meson V was chosen as in [24], Masses, total widths and 
branching ratios of the resonances were taken from [25]. The following value 
of the cross section at the 0 meson peak was obtained from the fit: 


cr(0 -> K° l K° s ) = (1376 ± 6 ± 23) nb, 
x 2 /n.d.f = 94.64/56 = 1.69. 

The value of the peak cross section is exactly equal to that from our 0 —» K^Kg 
study [22]. The relatively large value of the y 2 arises from the energy range 
above 1.13 GeV where, as seen from Fig. 2, most of the experimental points 
lie above the approximation curve. One of the possible explanations for the 
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Table 1 

The c.m. energy, number of events, integrated luminosity, detection efficiency, ra¬ 
diative correction, “dressed” cross section, vacuum polarization factor and “bare” 
cross section of the process e + e _ —> K®Kg. The first error in the cross section value 
is statistical and the second is systematic one. 


y/s, GeV 

N 

C, nb 1 

£ 

(1 + 5 ra( i ) 

a, nb 

|i - n (s )| 2 

cr°, nb 

1.050 

310.9T17.8 

117.7±0.7 

0.151 

0.996 

17.56=bl.00=b 1.97 

0.948 

16.65±0.95±1.87 

1.060 

124.1T11.3 

76.8±0.6 

0.146 

0.950 

11.65=1=1.06± 0.64 

0.953 

11.10=1=0.97=1=0.61 

1.070 

76.14=8.9 

81.8=10.5 

0.141 

0.926 

7.13=t0.84=t 0.35 

0.956 

6.82=t0.80±0.33 

1.080 

39.3T6.6 

59.6=t0.5 

0.136 

0.914 

5.30=t0.90=t 0.26 

0.958 

5.08=1=0.86=1=0.25 

1.090 

53.0T7.5 

85.2=tl.3 

0.132 

0.909 

5.18=t0.73=t 0.25 

0.959 

4.97=1=0.70=1=0.24 

1.100 

29.0T5.5 

57.5=t0.5 

0.128 

0.906 

4.35=t0.81=t 0.21 

0.961 

4.19=t0.78=b0.20 

1.110 

33.5T6.1 

83.3=t0.5 

0.124 

0.905 

3.58=t0.65=t 0.29 

0.962 

3.44±0.62±0.28 

1.120 

18.2^2 

59.4=t0.7 

0.120 

0.905 

2.82±°;!!±0.23 

0.963 

2.7li!;!!±0.22 

1.130 

14 7+ 4 - 2 

62.8=t0.5 

0.116 

0.905 

2.23t!;!!=b0.18 

0.963 

2.15i°;!!±0.18 

1.140 

25-51 li 

87.4=t0.5 

0.113 

0.905 

2.85^!;!!=b0.23 

0.964 

2.75i!;!!±0.22 

1.150 

17.8lt? 

63.2=t0.5 

0.109 

0.906 

2.85l!;^=t0.23 

0.965 

2.75i!; 6 7 !±0.21 

1.160 

17.2±£? 

113.0=t0.7 

0.106 

0.906 

1.58+!;^±0.14 

0.965 

i-52i!;i±o.i3 

1.170 

13-9111 

71.8=t0.5 

0.103 

0.907 

2.07±!;!!±0.17 

0.966 

2.00t!;!!±0.16 

1.180 

15.l!£g 

114.9=t0.7 

0.100 

0.907 

l-45l!;|±0.12 

0.966 

1-4018;£±0.12 

1.190 

20.91!;! 

132.5=tl.l 

0.097 

0.908 

1.79+!;^±0.15 

0.966 

1.72l°;i±0.14 

1.205 

19.71!;! 

270.1=tl.l 

0.094 

0.909 

0.85^!|!=b0.07 

0.967 

0.82t!;|!±0.07 

1.225 

26.91!;! 

286.7=tl.5 

0.089 

0.910 

1.16l!|!±0.11 

0.967 

1 . 12 +!;!!± 0.11 

1.251 

21-31!;! 

492.3=t2.3 

0.083 

0.912 

0.57+!;^±0.05 

0.968 

0.55t!;}!±0.05 

1.275 

18.71!;! 

468.6=11.5 

0.079 

0.913 

0.55^!;!!±0.05 

0.968 

0.53t!;l!±0.05 

1.296 

20 . 21 !;! 

589.2=t2.7 

0.075 

0.915 

0.50!:!;!!=b0.05 

0.968 

0.48^!;}!±0.05 

1.325 

11 . 21 !;! 

994.9=t3.1 

0.071 

0.916 

0.17+!;!!=b0.02 

0.968 

0.16+!;!!±0.02 

1.368 

20.71!;! 

1338.7=t3.4 

0.066 

0.919 

0.25^!;!!±0.02 

0.968 

CM 

O 

o 

-H 

00 t> 

oo 

o o 
+ 1 

CM 
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observed excess could be higher resonances contributing to the production of 
pairs of neutral kaons. This assumption is confirmed by the combined analysis 
of our results with those from the DM1 detector obtained in the energy range 
1.4 - 2.18 GeV [1], see Fig. 3. To describe the energy dependence of the cross 
section at y/s ~ 1.6 GeV, the following amplitude was added to the amplitudes 
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Fig. 2. The cross section of the process e + e~ —> K®Kg in the energy range 
\/s = 1.05 - 1.38 GeV, measured in different experiments. The curve is the Vec¬ 
tor Dominance Model prediction with the contributions of the p( 770), u>( 782) and 
0 ( 1020 ) mesons. 

of the p( 770), cu(782) and 0(1020) mesons: 


A x = 


m 7 x r 2 x a(X^K 0 L K° s )q3(ml) 


q 3 (m' x ) 


m 


x 


S + lyfsTx 


■ e l5x 


The parameters cr(X —> K®Kg), Mx, ^x,^x as well as the cross section at 
the 0 meson peak were obtained from the approximation: 


cr(0 -»• K° l K° s ) = (1375 ± 6 ± 23) nb, 
a(X -x K° l K ° s ) = 0.73 ± 0.33 nb, 

M x = 1623 ± 20 MeV/c 2 , 
T x = 139 ± 60 MeV, 

5 X = 160° ±42°, 
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Fig. 3. The cross section of the process e + e~ —► K®Kg in the energy range 
v/i = 1.05 - 2.2 GeV compared to the Vector Dominance Model predictions: p{ 770), 
uj( 782), 0(1020) mesons (the dotted curve) and p( 770), u( 782), 0(1020) and X (the 
solid curve). 


x 2 /n.d.f = 56.42/62 = 0.91. 

The value of the cross section at the 0 meson peak changes from 1376 nb to 
1375 nb after the contribution of a higher resonance is taken into account. 
This change agrees with the estimation of the systematic uncertainty of the 
cross section (of about 1.5 nb) due to the model dependence of the cross 
section value [26]. The approximation is shown in Fig. 3 by the solid line. 
The theoretical curve describing the contributions of the p(770), cu(782) and 
0(1020) mesons only, is shown in Fig. 3 by the dotted line. One can see that 
after adding the amplitude A x , the quality of the approximation of the exper¬ 
imental data is much better. The values obtained for the mass and width of 
the X state are consistent with those of the 0(1680) meson [25]. However, as 
noted above, both isovector and isoscalar states could contribute to the cross 
section of the neutral kaon production in the energy range 1.0 - 2.0 GeV. To 
identify unambiguously the nature of the observed enhancement as well as 
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Table 2 

Main sources of the systematic errors 


Source 

Contribution,% 

Selection criteria 

3.5-6 

Background subtraction 

2-8 

Luminosity 

2 

Detection efficiency 

2 

Radiative corrections 

1 

Total 

5-10 


to determine the relative weights of the 1=0 and 1=1 final states, one should 
simultaneously study both K®Kg and K + K~ final states with a significantly 
higher data sample, particularly in the energy range above 1.4 GeV. The de¬ 
tailed investigation of other final states of e + e~ annihilation in this energy 
range will be also needed to shed light on the spectroscopy of the light quark 
resonances in the vector sector. Such studies will be possible at the VEPP- 
2000 collider [27] now under construction at the Budker Institute of Nuclear 
Physics. 

The main sources of the systematic uncertainties are listed in Table 2. The 
uncertainty caused by the selection criteria was estimated by varying the cuts 
for pion momenta, acollinearity angle Acp and Z-coordinate of the vertex by 
one standard deviation. It showed that the cross section changed by 3.5% be¬ 
low 1.1 GeV and by 6 % above this energy. To estimate the uncertainty due 
to the background shape, the invariant mass distribution of two tracks was 
approximated in a narrow mass range from 400 to 600 MeV/c 2 , the parame¬ 
ters of the function describing background were determined and the obtained 
number of signal events was compared to that after the standard approxima¬ 
tion in the invariant mass range from 340 to 700 MeV/c 2 . The variation of 
the cross section was 2 % in the energy range 1.05-1.09 GeV smoothly rising 
to 8 % in the higher energy range. The systematic error of the luminosity is 
caused by the radiative corrections to the Bhabha scattering cross section as 
well as selection of Bhabha events [17]. The uncertainty of the detection ef¬ 
ficiency is dominated by the systematic error in the reconstruction efficiency 
and was estimated from the statistical error of the approximating constant 
in experiment [19]. The uncertainty on the radiative correction comes from 
the accuracy of the theoretical formulae used in the calculation ( 1 %) as well 
as from the missing momentum resolution. To take into account the influ¬ 
ence of the latter effect, selection criteria using this parameter were varied by 
one standard deviation. The resulting change was significant at 1.05 and 1.06 
GeV only where it was 10% and 2.5% respectively whereas at all other ener¬ 
gies it was negligible. The overall systematic uncertainty of the cross section is 
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obtained by adding the individual contributions in quadrature. It grows with 
energy from 5% in the energy range yfs = 1.05-1.09 GeV to 10% in the energy 
range y/s = 1.27-1.38 GeV. 


4 Conclusion 


Using 948±33 reconstructed events detected by the CMD-2, the cross sec¬ 
tion of the process e + e~ —> K Q L K° S was determined in the energy range from 
1.05 to 1.38 GeV. This is the most precise measurement of this cross section 
by now. It is shown that in the energy range yfs > 1.13 GeV the obtained 
energy dependence of the cross section could not be explained by the Vector 
Dominance Model with the contributions of the p(770), u;(782) and 0(1020) 
mesons only, so that higher resonances should be taken into account. The 
addition of an amplitude corresponding to a resonance with mass and width 
close to those of the 0(1680) meson and interfering with the p(770), u;(782) 
and 0(1020) mesons substantially improves the description of the observed en¬ 
ergy dependence. To obtain the detailed information about the spectroscopy 
of higher resonances, new experiments are needed which will be performed at 
the VEPP-2000 collider. 
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